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Abstract. We present local velocity measurements in emulsions under shear using heterodyne Dynamic 
Light Scattering. Two emulsions are studied: a dilute system of volume fraction <j> = 20 % and a concen- 
trated system with <f> — 75 %. Velocity profiles in both systems clearly show the presence of wall slip. We 
investigate the evolution of slip velocities as a function of shear stress and discuss the validity of the correc- 
tions for wall slip classically used in rheology. Focussing on the bulk flow, we show that the dilute system 
is Newtonian and that the concentrated emulsion is shear-thinning. In the latter case, the curvature of 
the velocity profiles is compatible with a shear-thinning exponent of 0.4 consistent with global rheological 
data. However, even if individual profiles can be accounted for by a power-law fluid (with or without a yield 
stress), we could not find a fixed set of parameters that would fit the whole range of applied shear rates. 
Our data thus raise the question of the definition of a global flow curve for such a concentrated system. 
These results show that local measurements are a crucial complement to standard rheological tools. They 
are discussed in light of recent works on soft glassy materials. 



PACS. 83.80.Iz Emulsions and foams - 83.50.Rp Wall slip and apparent slip 
yield stress - 83.85.Ei Optical methods; rheo-optics 
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1 Introduction 



Emulsions are mixtures of two immiscible fluids. They 
are usually formed by mechanical mixing and consist of 
droplets of one fluid dispersed in the other. They are 
metastable so that the mean size of the droplets tends 
to increase with time. Emulsions belong to a wide class 
■ of non-equilibrium systems such as foams that rearrange 
' and coarsen with time JlJ . The use of surfactants or poly- 
mers increases the characteristic time for coarsening of 
emulsions from a few seconds to a few years and make 
them very attractive for applications such as road build- 
ing, pharmacology or food processing. They constitute 
very interesting tools to entrap hazardous solvents, to 
carry additives or to change the rheological properties of 
a sample. Indeed, their rheological characteristics may be 
tuned by increasing the liquid fraction of the dispersed 
phase. 

When the liquid fraction is low, the droplets of the 
emulsion are Brownian and the emulsion is a Newtonian 
fluid whose viscosity is close to that of the continuous 
phase. When the liquid fraction is increased, less space is 
available to each particle. The droplets are pressed against 
each other, their shape are distorted, and energy is stored 
at their interfaces. For small deformations, the droplets 
resist shear elastically. Such elasticity results from the en- 
ergy stored by additional deformation of their shape in- 



duced by the applied strain. For large deformations, they 
flow irreversibly. 

Mason et al. || have shown that this change of regime 
occurs above a critical stress and that, for concentrated 
emulsions in the neighbourhood of this critical stress, the 
flow may become inhomogeneous. In Rcf. the homo- 
geneity of the flow was determined by painting a stripe 
on the exposed surface of the emulsion before shearing 
and by following its evolution under shear. Moreover, the 
analysis of the rheological data may be complicated by the 
existence of slippage of the emulsion at the walls . 

Due to all these difficulties, local measurements of the 
velocity arc needed to characterize the flow of the emulsion 
once it has yielded. Recently, Coussot et al. Q measured 
velocity profiles of various colloidal systems such as gels, 
laponitc solutions or industrial emulsions by imaging the 
flow held with Magnetic Resonance Imaging (MRI). They 
described the yield stress phenomenon as an abrupt tran- 
sition between two distinct states: a liquid-like state that 
flows, and a solid-like state that remains jammed and does 
not flow. The rheological behaviour of the liquid state was 
modelled by that of a power-law fluid. 

In this article, we present measurements of velocity 
profiles in sheared emulsions using Dynamic Light Scatter- 
ing (DLS) in the heterodyne geometry. Our setup allows us 
to perform both global rheological measurements and local 
velocity measurements simultaneously. We study two dif- 
ferent oil- in- water emulsions: a dilute emulsion of volume 
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fraction <j> = 20 % and a concentrated one at </> = 75 %. 
In both cases, we get a quantitative measure of wall slip 
and show the importance of this phenomenon. Wc point 
out that corrections for slippage like those proposed by 
Mooney and classically used in rheological experiments arc 
indeed well suited for concentrated emulsions but may fail 
for dilute systems. In the case of the concentrated emul- 
sion, we show the existence of two regimes separated by 
an abrupt transition between a jammed state and a flow- 
ing phase like Coussot et al However, in our case, no 
clear analytical flow curve could be obtained for the whole 
range of shear rates investigated in the present study. 

In the next Section of this article, we describe the 
method used for preparing the samples. The third Sec- 
tion briefly recalls our experimental technique and setup 
for measuring the local velocity of a complex fluid in a 
Couette flow Q . We then present the velocity profiles ob- 
tained in the two emulsions under shear. Finally, we dis- 
cuss both the slippage and the yield stress phenomena in 
our emulsions. 



2 Preparation of the emulsions under study 

We prepare a crude emulsion of fixed composition by gen- 
tly shear-mixing 400 g of silicone oil (Polydymethyl Silox- 
ane of viscosity 135 Pa.s from Rhodia), 30.3 g of glyc- 
erol, 30.3 g of water and 9.9 g of surfactant (Tetradecyl 
Trimethyl Ammonium Bromide from Aldrich). The aque- 
ous phase is thus composed of a 14 % wt. solution of TTAB 
in a 1:1 mixture of water and glycerol. The resulting poly- 
disperse viscoelastic premixed emulsion of large droplets 
is sheared within a narrow gap in a laminar regime |6[]. 




Fig. 1. Photograph of the emulsion under study. The con- 
centrated emulsion was diluted in order to visualize individual 
droplets using a microscope (magnification x50). The black 
bar represents 20 /xm. 

The mother emulsion obtained after shearing is com- 
posed of monodisperse silicone oil droplets of diameter 
2 fim in the water-glycerol-surfactant mixture (see Fig. [l]). 
Then, by dilution, we set the surfactant concentration 
within the water-glycerol phase to 1 % in mass and con- 
centrate the droplets by centrifugation to a volume frac- 
tion of 4> = 75 %, such that droplets are compressed one 
against another leading to flat films at each contact. We 



chose this surfactant concentration in order to ensure a 
good stability to the emulsion. This concentration is high 
enough to prevent coalescence and small enough to avoid 
flocculation by depletion. We prepare another emulsion 
with a volume fraction cf> = 20 % by diluting the concen- 
trated one in the water-glycerol-surfactant solution. The 
surfactant concentration is kept equal to 1 %. 

The use of a water-glycerol mixture allows us to match 
the optical index n of the continuous phase of the emul- 
sion to that of the silicone oil (n = 1.40). The obtained 
emulsions are thus nearly transparent and do not present 
any significant multiple light scattering. The emulsion is 
then carefully loaded into the transparent Couette cell of 
a rheometer (see below). For the whole range of applied 
shear rates, we check by optical microscopy that the struc- 
ture of the emulsion and the size of the droplets are not 
affected by shear and that no noticeable coalescence takes 
place. 

3 Experimental setup for heterodyne DLS in 
Couette flows 

In order to measure the local velocity of emulsions in Cou- 
ette flows, we used the experimental technique based on 
heterodyne DLS described at length in a related work H . 
In this Section, we briefly recall the basics of heterodyne 
DLS and the setup developed in our laboratory to per- 
form local velocity measurements. The reader is referred 
to Ref. § for more details on the optics of this setup and 
for a complete discussion on the resolution of heterodyne 
DLS. 




Fig. 2. Simplified heterodyne DLS setup. BS denotes a beam 
splitter, SF spatial filters, and C the device coupling optical 
fibers and used to perform the interference between the scat- 
tered light and the reference beam. 



Local velocity measurements using heterodyne DLS re- 
lie on the detection of the Dopplcr frequency shift asso- 
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ciated with the motion of the scatterers inside a small 
scattering volume V In classical heterodyne se- 

tups, light scattered by the sample under study is collected 
along a direction 9. The scattering volume is defined as the 
intersection between the incident beam and the scattered 
beam. The scattered electric field is then made to inter- 
fere with a reference beam. Finally, light resulting from the 
interference is sent to a photomultiplicr tube (PMT) and 
the auto-correlation function of the intensity is computed 
using an electronic correlator. The originality of our setup 
lies in the use of single mode fibers to collect scattered 
light, to perform the interference, and to carry light to 
the PMT. This allows for great flexibility when choosing 
the scattering angle. Figure || sketches the main features 
of such a setup. 

When the scattering volume V is submitted to a shear 
flow, it can be shown that the correlation function C(t) 
is an oscillating function of the time lag t modulated by 
a slowly decreasing envelope. The frequency of the oscil- 
lations in C(t) is exactly the Doppler shift q ■ v, where 
q is the scattering wavevector and v is the local velocity 
averaged over the size of the scattering volume which is 
about 100 /im in our experiments. Figure || shows a typi- 
cal correlation function measured on a sheared emulsion. 
The frequency shift q • v is recovered by interpolating a 
portion of C(r) and looking for the zero crossings. Error 
bars on such measurements are obtained by varying the 
number of zeros taken into account in the analysis. Typical 
uncertainties arc about 5 %. 




Fig. 3. Experimental heterodyne correlation function (•) 
recorded on the concentrated emulsion at 7 = 10 s . The 
thick line shows the portion of the heterodyne function from 
which the frequency shift q ■ v is calculated. 



Couette flows are generated between two cylinders in a 
transparent Mooncy-Couette cell. The radius of the rotat- 
ing inner cylinder ( "rotor" ) is R t = 22 mm and the radius 
of the fixed outer cylinder ( "stator" ) is R 2 = 25 mm leav- 
ing a gap e = R 2 — Ri = 3 mm between the two cylinders. 
The length of the cylinders is H — 30 mm. Both cylinders 
are in Plexiglas and have smooth surfaces so that wall 
slip is likely to occur in our experiments. The rotation of 



the moving cylinder is imposed by a classical rheometer 
(TA Instruments AR1000) operated in the "imposed shear 
rate" mode. This rheometer measures both the torque r 
on the inner cylinder and its angular velocity fi in real- 
time. From fl and -T, it also computes a global shear rate 
< 7 > and shear stress < a > using the formulas: 



< 7 > 



R, + R., 



< a > = 



Rl- 
R 2 



- Ri 



AnHR 2 R 2 



r. 



(i) 

(2) 



In the following, the results will be presented in terms of 
< 7 > and < a > which will simply be noted 7 and er, as 
is usually done in rheology. Using Eq. (|l|), the velocity of 
the moving wall v = Rift is given by: 



Ri{Ri + R2) 

R? + R 2 



(■■1) 



Thus, in the Couette geometry used in the present work, 
one has to keep in mind that v a ps 0.93 je so that v differs 
from je by about 7 %. 

The temperature of the sample is controlled within 
±0.1°C using a water circulation around the cell. The 
rheometer sits on a mechanical table whose displacements 
are controlled by a computer. Three mechanical actuators 
allow us to move the rheometer in the x, y, and z directions 
with a precision of 1 fim. Once y and z are set so that 
the incident beam is normal to the cell surface, velocity 
profiles are measured by moving the mechanical table in 
the x direction by steps of 50 or 100 fj,m. 

As discussed in Ref. || , going from q- v as a function of 
the table position xt to the velocity profile v(x) (where x 
is the radial position of the scattering volume V within the 
sample) requires a careful calibration procedure. Indeed, 
refraction effects due to the curved geometry of the Cou- 
ette cell lead to significant differences between the angle 9i 
imposed by the operator and the actual scattering angle 9 
in the sample, and between the table displacement 5x and 
the actual displacement 6x' of V. By using a Newtonian 
suspension of latex spheres in a water-glycerol mixture 
(optical index n = 1.40) at various known shear rates, we 
showed that for 9i = 35° , 9 = fg9i and Sx' = f x Sx with 
fe « 0.79 and f x « 1.13. 

When loading the cell with the emulsions, great care is 
taken to ensure that the position of the Couette cell docs 
not change between the calibration with the latex suspen- 
sion and the actual experiments on emulsions. Since the 
optical index of the fluid is n — 1.40 in both cases, the 
above values of f x and fg found with the latex suspension 
are used to convert raw data (q • v)(xt) into velocity pro- 
files v(x) in the case of our emulsions. Finally, we chose to 
accumulate the correlation functions over 1 min so that 
a full velocity profile with a resolution of 100 /im takes 
about 30 min to complete. Note that our setup allows us 
to record rheological data and the local velocity at the 
same time on the same sample. In all the cases described 
below, we checked that steady state was reached from the 
curves ^(t) and cr(t). 
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4 Velocity profiles in sheared emulsions 
4.1 Dilute emulsion ((f) = 20 %) 




x (mm) x (mm) 

Fig. 4. Velocity profiles in the 20 % emulsion for (a) 7 = 
0.1 s -1 and (b) 7 = 0.2 s _1 . The solid curves are the best 
fits to the data using Eq. with n = 1. The dashed lines 
represent the velocity profiles expected for a Newtonian fluid. 




x (mm) x (mm) 

Fig. 5. Velocity profiles in the 20 % emulsion for (a) 7 = 
0.4 s _1 and (b) 7 = 1 s _1 . The solid curves are the best fits to 
the data using Eq. ([[5]) with n = 1. The dashed lines represent 
the velocity profiles expected for a Newtonian fluid. 

Figures ^-|^ present the velocity profiles measured in 
the dilute emulsion of liquid fraction <f> = 20 % for various 
imposed shear rates ranging from 0.1 s _1 to 15 s" 1 . The 
velocity profiles expected for a Newtonian fluid are shown 
as dashed lines. The velocities v(x) have been rescaled by 
the inner wall velocity v . In these figures, the upper left 
corner (at x — 0) thus always corresponds to the velocity 
of the inner wall moving at v , and the lower right corner 
(at x = e = 3 mm) to the fixed outer wall (v = 0). First, 
one can notice that the experimental velocity is always 
different from zero at the stator and different from the 
rotor velocity at the rotor. This means that significant 
slippage occurs and remains of the order of 10 % ofv . We 
define the slip velocity at the rotor v Sl = v — v t as the 
difference between the rotor velocity v and the emulsion 
velocity v 1 at the rotor and the slip velocity v S2 = v 2 at 
the stator as the emulsion velocity v 2 at the stator. 

Moreover, all the velocity profiles measured in the di- 
lute emulsion are linear. To better illustrate this point, 
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x (mm) x (mm) 

Fig. 6. Velocity profiles in the 20 % emulsion for (a) 7 = 2 s _1 
and (b) 7 = 5 s _1 . The solid curves are the best fits to the 
data using Eq. ( |l5| ) with n = 1. The dashed lines represent the 
velocity profiles expected for a Newtonian fluid. 




x (mm) x (mm) 

Fig. 7. Velocity profiles in the 20 % emulsion for (a) 7 = 10 s _1 
and (b) 7 = 15 s . The solid curves are the best fits to the 
data using Eq. ( |l5| ) with n = 1. The dashed lines represent the 
velocity profiles expected for a Newtonian fluid. 




x (mm) 



Fig. 8. Normalized velocity profiles measured in the 20 % 
emulsion for 7 = 0.1 s" 1 (V), 0.2 s" 1 (A), 0.4 s" 1 (+), 
1 s" 1 (x), 2 s" 1 (o), 5 s" 1 (□), 10 s" 1 (•), 15 s" 1 (o), 
and 20 s _1 (*). Each profile is normalized according to v n = 
(v — v 2 )/(v 1 — v 2 ). The dashed line is the result expected for 
a Newtonian fluid. Inset: the average of the previous normal- 
ized profiles < v n > (•) compared to the Newtonian prediction 
(dashed line). 
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Fig. ^| presents the normalized velocity v n for various ap- 
plied shear rates. We define v n experimentally by: 



v n (x) : 



v (x) - V2 
vi - v 2 



(4) 



All normalized data coincide and the average master curve 
< v n > is almost perfectly linear across the cell gap (see 
inset of Fig. ||). These qualitative results show that, al- 
though subject to small but measurable wall slip, the di- 
lute emulsion flows like a Newtonian fluid in the range of 
shear rates under study. 



4.2 Concentrated emulsion (</> = 75 %) 





Fig. 11. Velocity profiles in the 75 % emulsion for (a) 7 = 2 s _1 
and (b) 7 = 5 s _1 . The solid curves are the best fits to the data 
using Eq. ( |l5| ) with n = 0.4. The dashed lines represent the 
velocity profiles expected for a Newtonian fluid. The dotted 
lines correspond to the case of a Newtonian fluid that would 
slip at the cell walls with the same velocities as the emulsion. 
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Fig. 9. Velocity profiles in the 75 % emulsion for (a) 7 = 
O.f s _1 and (b) 7 = 0.2 s . The solid lines are linear fits 
to the data. The dashed curves represent the velocity profiles 
expected for a Newtonian fluid. 
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Fig. 10. Velocity profiles in the 75 % emulsion for (a) 7 = 
0.4 s _1 and (b) 7 = 1 s - . The solid curves are the best fits 
to the data using Eq. ji5| ) with n = 0.4. The dashed lines 
represent the velocity profiles expected for a Newtonian fluid. 
The dotted lines correspond to the case of a Newtonian fluid 
that would slip at the cell walls with the same velocities as the 
emulsion. 

Figures p|-|l2] present the velocity profiles in the con- 
centrated emulsion of liquid fraction <f> = 75 % for the 
previous range of imposed shear rates. Two kinds of be- 
haviours are encountered. In the case of low applied shear 





Fig. 12. Velocity profiles in the 75 % emulsion for (a) 7 = 
10 s~ and (b) 7 = 15 s" 1 . The solid curves are the best fits 
to the data using Eq. (^) with n = 0.4. The dashed lines 
represent the velocity profiles expected for a Newtonian fluid. 
The dotted lines correspond to the case of a Newtonian fluid 
that would slip at the cell walls with the same velocities as the 
emulsion. 



rates (0.1 and 0.2 s , see Fig. |9|), the emulsion moves 
around in the gap as a solid body with angular velocity 
uj = 0.042 rad.s" 1 for 7 = 0.1 s" 1 and u = 0.084 rad.s" 1 
for 7 = 0.2 s _1 . In this regime, all the material actually 
experiencing shear is located inside thin boundary layers 
at the two walls. These boundary layers cannot be re- 
solved with our spatial resolution of about 100 fim. When 
the applied shear rate is increased, the velocity profiles 
change. Wall slip remains very important but the emul- 
sion no longer behaves as a solid body. Shear occurs not 
only in the boundary layers but also in the bulk material. 

The velocity profiles arc no longer linear but rather 
curved (see Fig. O for instance). Moreover, as shown in 
Fig. [e| the data for the normalized velocity v n (x) still 
collapse rather well on a single curve provided the applied 
shear rate is greater than 0.2 s _1 . The curvature of the 
mean normalized profile < v n > is clearly visible in the 
inset of Fig. [13[ This result means that the non-Newtonian 
features of our concentrated emulsion are strong enough 
to show up at a local level even in a gap as small as 3 mm. 
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Fig. 13. Normalized velocity profiles measured in the 75 % 
emulsion for j = 0.4 s" 1 (+), 1 s" 1 (x), 2 s _1 (o), 5 s _1 (□), 
10 s _1 (•), and 15 s _1 (o). Each profile is normalized accord- 
ing to v n = (v — v 2 )/(vi — v 2 ). The dashed line is the result 
expected for a Newtonian fluid. Inset: the average of the previ- 
ous normalized profiles < v n > (•) compared to the Newtonian 
prediction (dashed line). 

In the next Section, we first discuss the evolution of the 
slip velocities and then analyze the features of individual 
velocity profiles in more depth. 

5 Local rheology of emulsions in a Couette 
flow 



Fig. 14. Wall slip velocities for the 75 % emulsion (•) and the 
20 % emulsion (o) vs. shear rate, (a) Normalized slip velocity 
v si /v at the moving inner wall (x = 0). The solid line is the 
power law 0.35 -y -0 ' 4 . The dotted line shows the mean value 
v si = 0.048 v for the 20 % emulsion, (b) Normalized slip ve- 
locity v S2 /vo at the fixed outer wall (x = 3 mm). The solid line 
is the power law 0.20 7 -0 ' 4 . The dotted line shows the mean 
value v S2 — 0.10 Wo for the 20 % emulsion. 

In the dilute system, the normalized slip velocities at 
the rotor v si /v and at the stator v S2 /v are independent 
of the applied shear rate as can be checked on Fig. [l4| 
This means that v Si ~ 7 and since the dilute emulsion has 
a Newtonian behaviour, one expects v Si ~ (Tj. This scaling 
law is indeed observed on Fig. [15|(a). 



5.1 Wall-slip and wetting films 

For both the dilute and the concentrated emulsions, slip- 
page clearly occurs. The simplest picture one could imag- 
ine to describe this situation is the existence of two very 
thin films of continuous phase near the walls and a bulk 
material with the original concentration in between jlO 1 1 . 
From our measurements, it is possible to access the thick- 
ness h of the thin film where shear occurs. Indeed, assum- 
ing a slip velocity v s within a thin liquid film of viscosity 
rjf, we may write: a = r)fj s = rjfv s /h, where 7 S = v s /h is 
the shear rate inside the film. Using the subscript i = 1 or 
2 to denote the film at the rotor or at the stator respec- 
tively, 

hi = for i = 1 or 2 , (5) 

where <Ji is the shear stress near wall number i. Note that, 
in the Couette geometry, the values of Oi at the walls are 
linked to the global stress a indicated by the rheometer 
according to: 

ai = WTm^ (6) 

where J = 2 (resp. j = 1) when i = 1 (rcsp. i = 2). For 
each value of 7, the steady-state value of a imposed by the 
rheometer is recorded and <j 1 and a 2 are computed from 
Eq. (||). In our case, this leads to significant deviations 
from a: a 1 = 1.13 a and a 2 = 0.87 a. 




0.3 0.6 0.9 200 400 600 

a. (Pa) a. (Pa) 



Fig. 15. Slip velocities vs. shear stress: v si vs. ai at the rotor 
(•) and v S2 vs. a 2 at the stator (o). (a) For the 20 % emulsion. 
The straight lines correspond to Vi = 33<7i and to v 2 = 8.0a 2 . 
(b) For the 75 % emulsion. The solid line is Vi = 1.5.10" 5 of. 

Assuming that the films contain only water and glyc- 
erol, one has 77/ = 0.01 Pa.s. The thicknesses of the films 
h 1 and h 2 calculated from Eq. (H) are shown as functions 
of <j 1 and <t 2 in Fig. fl6](a). One gets h 1 w 30 /xm at the 
rotor and h 2 w 80 (im at the stator independent of the 
shear stress, which is again consistent with a Newtonian 
behaviour. Remarkably, these values are much higher than 
the diameter of the droplets (2 /xm), meaning that deple- 
tion in the dilute emulsion extends very far from the walls. 
Various explanations may be proposed for such a large de- 
pletion effect. 
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a. (Pa) a. (Pa) 

Fig. 16. Thickness of the sliding layers vs. shear stress: hi 
vs. <7i at the rotor (•) and h 2 vs. <r 2 at the stator (o). (a) 
For the 20 % emulsion. The dashed lines are h ± = 34 /xm 
and h 2 — 83 fjxa. (b) For the 75 % emulsion. The solid line 
corresponds to hi = 1.5.10 - Uj. 

First, depletion on such long distances might be in- 
duced by cumulative effects of inertial effect and of cross- 
stream migration. Indeed, Goldsmith and Mason reported 
rapid migration from walls of neutrally buoyant and de- 
formable droplets in pipe and Couette flows |l2|. The au- 
thors suggest that such a migration could be initiated by 
the presence of the wall and of non-uniform velocity gra- 
dients. A net force acting to push the drop to a lower 
gradient is created. 

Second, it is important to note that the viscosity of 
the emulsion varies dramatically with oil concentration. 
For example, the viscosity of the 20 % emulsion indicated 
by the rheometer is 0.037 Pa.s. When diluted to a volume 
fraction <j) = 10 %, the viscosity decreases to 0.017 Pa.s. 
Finally, at <j) = 5 %, the viscosity is equal to 0.011 Pa.s, 
very close to that of the continuous phase. This means that 
the lubricating layers could also be viewed as emulsion 
films with a concentration gradient, for example from </> = 
20 % in the bulk to <fi < 5 % near the walls, rather than a 
film of pure continuous phase. This would reduce greatly 
the magnitude of the depletion and could extend its effect 
over large distances. 

Finally, the thicknesses of the films are different at the 
rotor and at the stator. Slippage occurs preferentially at 
the stator (about 10 % of v ) and remains about two times 
smaller at the moving inner wall (about 5 % of v ). Again, 
since the oil density is slightly below that of the water- 
glycerol mixture, this might be due to inertial effects that 
lead to different compositions of the two sliding layers. 

5.1.2 Concentrated emulsion ((j> = 75 %) 

In the concentrated system, the thicknesses of the lubri- 
cating layers are much smaller because for a given shear 
rate, the stress in the emulsion is much larger. It is clear 
that cross-migration and inertial effects are completely 
suppressed in this case. As shown in Fig JT3(b) , the thick- 
ness of the films calculated from Eq. dSjincreases from 
h 1 10 nm at low shear stress to h x « 90 nm at high 
shear stress. These thicknesses are of the order of magni- 
tude of that of the thin liquid films of continuous phase 



that lie between the droplets. Such results are in qualita- 
tive agreement with previous measurements performed by 
Princcn using direct visualization of a stripe painted on 
the emulsion surface [flo|| . 

Several reasons may be given for the increase of h as a 
function of shear stress, (i) h may increase as a result of the 
cell walls dragging bulk continuous phase from the pockets 
between the drops inside the films, (ii) The viscosity of the 
wetting film may decrease due to an increased mobility of 
the drop/film interface as the stress increases. 

Like in the dilute emulsion, the phenomenon of wall 
slip may not seem symmetrical when one plots v Si /v vs. 
7 (see Fig. |lj). This time, for a given shear rate, the slip 
velocity is about 50 % larger at the rotor than at the sta- 
tor. However when plotted against the local shear stresses 
at the walls, the data v Sl (<Ti) and v S2 (a 2 ) collapse on each 
other as can be seen on Fig. ||(b). One can thus interpret 
the difference between v Sl (7) and v S2 (7) as a consequence 
of the difference between a 1 and tr 2 due to the inhomogene- 
ity of the shear stress in the gap. Indeed, in our Couette 
cell, the shear stresses at the rotor and at the stator differ 
by about 25 %. One then recovers a 50 % difference in the 
slip velocities provided v Si ~ of. This scaling behaviour 
is evidenced in Fig. |l5|(b) where a parabola is seen to fit 
nicely all the v s (a) data. A direct consequence of Eq. (||) 
is that h 1 (a 1 ) = h 2 (a 2 ) and that h(a) ~ a (see Fig. |l6|(b)). 

5.1.3 Discussion on usual corrections for wall slip 

The previous results are important because they provide 
a confirmation that the corrections for wall slip proposed 
by Yoshimura and Prud'homme p| or by Kiljafiski [ [l4| 
based on those introduced by Mooncy jL5| are relevant 
for concentrated emulsions. Indeed, these authors assumed 
(i) that the slip velocity and the thickness of sliding lay- 
ers are functions of the shear stress only and (ii) that 
the slip velocities at the rotor and at the stator are the 
same functions of a. Then, from rheological experiments 
performed in different geometries (by varying the height 
p3| or the gap Q of the Couette cell), corrections for 
wall slip are computed, that yield an indirect estimate of 
v s (a). The present results based on direct measurements 
of v s at the two walls show that v Sl (a 1 ) = v S2 {<r 2 ) and 
reveal a quadratic behaviour of v s vs. a consistent with 
Refs. HQ. 

However, the indirect corrections of Refs. fll3,|l4i arc 
clearly not valid for our dilute emulsion since FigT[15[snows 
that v Sl (a 1 ) 7^ v S2 (a 2 ). In this case, local measurements 
are required to estimate accurately the influence of wall 
slip on rheological data. 

5.2 Local rheology vs. global rheology 

5.2.1 Effective shear rate 

Using the slip velocities estimated from the velocity pro- 
files, it is easy to define an effective shear rate in the emul- 
sion by 7eff = (max(w)— min(u))/e = (v 1 —v 2 )/e. However, 
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in order to take into account the stress inhomogencity in 
the Couette cell and thus get a quantitative comparison 
between this effective shear rate and the global shear rate 
7 given by the rhcomctcr. we have to define 7 e g consis- 
tently with Eqs. (0) and (|) by: 



7off = 



R 1 + R 2 
R 1 (R 1 + R 2 



(7) 



In Fig. [I?], 7 ff is plotted against the imposed shear rate 
7 for the two emulsions. Surprisingly both data sets are 
very well fitted by straight lines with the exact same slope 
0.85 (provided 7 > 1 s _1 for the concentrated emulsion, 
see inset of Fig. |f*7| ). For the dilute Newtonian emulsion, 
7eff « 0.857 oniv confirms that the emulsion undergoes a 
constant wall slip of about 15 % independent of 7. This 
implies that due to wall slip, the value of the viscosity 
V = a H given by the rheometer is underestimated by 
about 15 % so that the actual viscosity is 77 ~ 0.043 Pa.s 
instead of 0.037 Pa.s. The fact that %s « 0.857 + cs t also 
for the 75 % emulsion seems to indicate that this value of 
0.85 is somehow characteristic of the oil droplets (radius, 
surface tension, etc.). Further studies should focus on this 
specific point. 




Shear rate (s ) 



Fig. 17. Effective shear rate in the 75 % emulsion (•) and 
in the 20 % emulsion (o) vs. imposed shear rate. The solid 
line is the best linear fit of the 75 % data for 7 > 1 s~ : 
7eff = 0.85 7 ~~ 0.61. The dotted line is the best linear fit of 
the 20 % data: 7 off = 0.85 7 + 0.005. Inset: blow-up of the low 
shear regime. 



The emulsion sample was the same as the one used for 
the velocity profile measurements. Striated surfaces are 
supposed to minimize wall slip whereas the measurements 
obtained with the smooth walls lead to the values of a and 
7 used up to now. 

A constant stress is applied on the sample for ten min- 
utes and the steady-state shear rate is recorded. The flow 
curves are obtained point by point from the lowest stress 
upwards. We checked that measurements from the highest 
stress downwards give the same results. The two steady- 
state flow curves (7(7) obtained for the concentrated emul- 
sion are compared on logarithmic scales in Fig. The 
two curves are clearly different . If one considers only the 
measurements obtained with the Couette geometry, it is 
not obvious whether the emulsion has a yield stress or not. 
Indeed, the yield stress inferred from this data set is very 
small: a k, 4 Pa. On the other hand, measurements ob- 
tained with the striated cone-and-plate geometry clearly 
present a yield stress of about 70 Pa. 




5.2.2 Global flow curve(s) 



Shear rate (s ) 



Fig. 18. Rheological flow curves a vs. 7 for the 75 % emulsion 
measured in a striated cone-and-plate geometry (T) and in the 
Plexiglas Couette cell with smooth walls and e = 3 mm (o). 
The solid line is a = 70 + 120 7 ' 4 and the equation of the 
dotted line is a = 4 + 220 7 4 . 



However, even when using striated surfaces, global rhe- 
ological measurements may still be affected by small slip- 
page or by fracture behaviour at shear rates close to the 
yield point j2|. The determination of a precise value for 
the yield stress is thus difficult with standard rheological 
tools. A way around the problem is to measure the lo- 
cal velocity and plot the shear stress as a function of the 
effective shear rate 7 e ff in the sample. 



In the particular case of the concentrated system, it is 
interesting to compare the classical rheological data 
to the "effective flow curve" a(%s) that can be estimated 
from the local velocity measurements. 

Classical steady-state flow measurements were carried 
out using both a cone-and-plate geometry with striated 
surfaces (cone angle 2° and radius 20 mm) and the Cou- 
ette cell with smooth walls used throughout this study. 



5.2.3 Local flow curve 

As seen in Fig. [l9| where the "effective flow curve" cr(7 ff) 
is compared to the global rheological data obtained in the 
Couette cell, a(%s) yields yet another flow curve, different 
from the curves of Fig. [f8|. Now, the two points where the 
emulsion does not flow lie on the axis 7 = and the data 
corresponding to 7 > 0.4 s" 1 seem to point to a higher 
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value of the yield stress around 90 Pa. As we shall sec 
below, the present data do not allow us to draw definite 
conclusions about the value of the yield stress nor even 
about the existence of a simple analytic flow curve for the 
concentrated emulsion. 




Shear rate (s ) 



"local flow curve" (7(2), c(x)) mainly arises from the ex- 
perimental estimation of the derivative in Eq. (@). For each 
shear rate, each cluster of points (j(x),a(x)) was fitted 
by a straight line. The resulting segments are plotted in 
logarithmic scales in Fig. [l^. It can be seen that (7 c ff , ci) 
always fall close to these first order approximations of the 
local flow curve. Note that due to the inhomogeneity of the 
stress in the Couette cell, the shear rates explored locally 
continuously cover the whole range 0.03 s _1 -20 s _1 . 

However, a jump in the local behaviour can be clearly 
identified between 7 = 1 s -1 and 7 = 5 s _1 and the local 
flow curve does not appear to be continuous at larger shear 
rates. Thus, within the precision of our data, it looks as 
though the local flow behaviour differs from one imposed 
shear rate to another, at least for 7 > Is. These results 
indicate that the very notion of a simple analytic flow 
curve may become questionable at high shear rates for our 
concentrated emulsion. In the next paragraph, we try to 
make this point clearer by focussing on the precise shape 
of the velocity profiles. 



Fig. 19. The "effective flow curve" a vs. 7 c f[ (•) and the "local 
flow curve" o(x) vs. 7(x) (dots) for the 75 % emulsion com- 
pared to the rheological data obtained in the Plexiglas Couette 
cell (o). The solid segments are linear fits of the local data o(x) 
vs. 7(2;). 



Since the velocity profiles are curved, one may argue 
that plotting the effective shear rate %g against the shear 
stress given by the rheometer is not such an accurate es- 
timate of the flow curve. However, using local measure- 
ments, one can go one step further and estimate the local 
shear rate 7(x) by differentiating the velocity profile v(x) 
according to: 

This was done directly on the raw velocity data with 
a first-order approximation for the derivative: dv/dx ss 
5v/Sx = (v(x + 8x) — v(x))/Sx. The minus sign in Eq. (||) 
accounts for the orientation of the x axis along decreasing 
velocities. 

Moreover, if the flow is stationary and the sample does 
not present any fracture, then Vc = and the local shear 
stress reads: 

^ = °'(^- < 9 > 

where R ± is the radius of the rotor and o x = <r{x = 0) 
the shear stress applied by the rheometer on the rotor. 
Note that this equation was already implicitely used at 
the walls when writing Eq. (||). For each velocity profile 
and for each value of x, the point (j(x). cr(x)) obtained by 
Eqs. (§) and (|) is plotted in Fig. |lg 

We check that the data (7(2;), o-(x)) are scattered around 
the effective flow curve cr(7 c ff). The dispersion seen in this 



5.3 Trying to account for the curvature of the velocity 
profiles 

5.3.1 Testing the flow behaviour a = A^ n 

Using MRI, Coussot et al. Q found that above some criti- 
cal stress <7 C , the velocity profiles of various flowing concen- 
trated materials were well described by a local power-law 
behaviour: 

7(2) = for a < g c , (10) 
a(x) = A^{x) n for a > a c . (11) 

Those two equations describe a flow behaviour with two 
separate branches. Equation ( [To| ) represents the jammed 
state and holds below the yield point a c . Equation ( [Tl| ) 
corresponds to a power-law fluid. The case n — 1 and 
<7 C = corresponds to the case of a Newtonian fluid (A 
is then the usual viscosity 77). n < 1 (resp. n > 1) means 
that the fluid is shear-thinning (resp. shear-thickening). 
This exponent n controls the curvature of the velocity 
profiles v(x). In light of the results of Ref. j|, wc wanted 
to test whether Eq. ([ll]) was able to account for the shape 
of our velocity profiles when the emulsion flows i.e. for 
7 > 0.2 s- 1 . 

Theoretical velocity profiles. Using Eqs. (||), @ and fj"l|), 

it is straightforward to calculate the flow profile above the 
yield point: 

R t +x 

v { x ) = v 2 — + 

«-i(5(t)r((^r->). 

(12) 
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and the velocities of the emulsion at the two walls are 
linked by: 

/ n \ l/n 




As before, a 1 denotes the shear stress at the moving in- 
ner wall: a 1 = a{x = 0). Using Eqs. (O), (O), and the 
definition of Eq. (0), the normalized flow profile is given 
by: 



\Ri+x ) 



2/n 



1 



. 2/n 

Ei) _i 
Ri ' 



(14) 



The previous approximation results from the assumption 
that v 1 /R 2 <C j{x) everywhere in the cell gap. Since v t < 
v o, i{x) < v /e, and R 2 — 10 e, Eq. ([IJ) holds with an 
accuracy that is always better than the experimental un- 
certainty of about 5 %. 



Experimental determination of the exponent n. Equa- 
tion ( |l4| ) shows that in the framework of Eq. ([ll]), fo- 
cussing on v n (x) allows us to remove the contributions 
of both the slip velocity v 2 and the prefactor A to the 
velocity profiles, so that the only free parameter is n. 
Figure 2(j(a) presents theoretical profiles obtained with 
Eq. (|l4) for R 1 = 22 mm, R 2 = 25 mm, and values of n 
ranging from 0.2 to 1. The inset shows that velocity pro- 
files can hardly be distinguished from one another when 
0.6 < n < 1. In Fig. ^(b), we plotted the master curve 
< Vn (x) > found for the concentrated emulsion (inset of 
Fig. [l3|) together with theoretical profiles for n = 0.3, 0.4, 
and 0.5 given by Eq. fll4|). The experimental data and 
their error bars fall between the theoretical profiles with 
ri = 0.3 and n = 0.5 so that we have n = 0.4 ± 0.1 in 
the case of the concentrated emulsion. This value of n is 
compatible with both the model of Princcn j[6| that pre- 
dicts an exponent n= 1/2 and that of Berthicr et al. JL7J 
which yields n = 1/3 close to the glass transition in soft 
materials. 

When performed on the dilute emulsion, the same anal- 
ysis yields n = 1 ± 0.1. In this case, the very small curva- 
ture of v n {x) due to the variation of a(x) across the gap 
cannot be detected and the Newtonian profile is undis- 
tinguishablc from a straight line, at least in a Couette 
cell with R 1 = 22 mm and e = 3 mm. On the other hand, 
when the fluid is shear-thinning enough, say when n < 0.6, 
the effect of the shear stress non-uniformity is enhanced 
and curvature becomes quantitatively measurable even in 
a Couette cell with a rather small gap. All the velocity 
profiles obtained with the concentrated emulsion are con- 
sistent with a single value of n = 0.4 ± 0.1 and more ex- 
periments, perhaps in wider gaps, could help discriminate 
between n =1/3 and n = 1/2. However, for the model 
to be fully applicable, one still has to check that Eq. ( pr| ) 
holds for a single value of the parameter A over the whole 
range of investigated shear rates. 




(b) 



1 2 
x (mm) 



Fig. 20. (a) Theoretical velocity profiles obtained with 
Eq. @ for n = 0.2 (•), 0.4 (□), 0.6 (x), 0.8 (o), and 1 (+). 
Inset: blow-up of the middle of the gap. (b) Average normal- 
ized profiles for the 75 % emulsion < v n > (•) and the non- 
Newtonian predictions of Eq. ( |l4| ) with n = 0.3 (lower dotted 
line), n = 0.4 (solide line), and n = 0.5 (upper dotted line). 
Inset: blow-up of the middle of the gap. 



Experimental determination of the parameter A and 
failure of the model. Once n is known, it is possible to es- 
timate the prefactor A from our experimental data. More 
precisely, the data in Figs. f|-(7|and |lo|-|l2l were fitted using 
a least-square algorithm by: 



v(x) = B 



Rx 



R* 



C(R ± 




(15) 

where B and C are free parameters and n was fixed to n = 
1 for the 20 % emulsion and n = 0.4 for the 75 % emulsion. 
The best fits are plotted in solid lines in the corresponding 
figures. In all cases, those fits describe the experimental 
data very closely. Note that for the concentrated emulsion 
at 7 — 0.1 s" 1 and 0.2 s _1 (Fig. ||), the data were fitted 
by straight lines since the emulsion undergoes solid-body 
rotation in this low shear regime. 

Moreover, if Eq. (|l^ ) holds, one expects: 



B = v 2 and C = 




(16) 



Hence, if the power-law behaviour of Eq. (|Tl|) is verified, 
one should have: 



R 
R, 



2C 
n 



1± 

A 



(17) 



Thus, an important self-consistency check of the model 
can be performed by plotting D against the shear stress at 
the rotor <j x . As shown in Fig. ^l|(a), the data for the dilute 
emulsion are very well fitted by D = a x /A with A ps 0.042. 
This is consistent with the Newtonian behaviour of the 
20 % emulsion and with the value found for its viscosity 
from rheological measurements when corrected for wall 
slip (r) = 0.043 Pa.s). 

As for the concentrated emulsion, the behaviour is 
quite different. Figure |l|(b) shows that the parameter D 
no longer behaves linearly with <r l and, more important, 
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0.3 0.6 0.9 200 400 600 



o, (Pa) Oj (Pa) 

Fig. 21. The parameter D vs. a (see text), (a) For the 20 % 
emulsion (o), the best linear fit yields A — 0.042. (b) For the 
75 % emulsion (•), the best linear fit for 150 < <Ti < 350 Pa 
(i.e. 0.4 < -y < 1 s _1 ) does not go to zero when a x — ► 0. 

that D no longer goes to zero when a 1 goes to zero. In- 
stead, D(a 1 ) intersects the axis D = at a shear stress a x 
of about 100 Pa which corresponds to a w 90 Pa (using 
Eq. (^)). This means that, when calculated by A = <r 1 /D 
according to Eq. (|l7|), A varies from 407 to 185 when 7 in- 
creases from 0.4 to 15 s . This strong variation of A with 
a shows that, even though the individual fits by Eq. ( |l5| ) 
are all very good, our data are not compatible with the 
flow behaviour given by Eq. ( |li"| ) with a single constant 
A. In particular, Eq. (|ll| ) may well hold locally and for a 
single value of n, but A has to be allowed to vary with a in 
order to account for the experimental data. Thus, in the 
case of our concentrated emulsion, the local model used 
in Ref. [|| fails to describe the whole range of investigated 
shear stresses with a single exponent n = 0.4 and a single 
value of A. 



shape when a 3> er , whereas the curvature of our experi- 
mental profiles does not vary noticeably (see Fig. [l3| ). 

We integrated Eq. (|l9| ) numerically and looked for a set 
of parameters v tl er , A, and n that would describe all the 
different profiles. Since v ± corresponds to slip at the rotor, 
the values v 1 (a) found earlier in this paper were forced 
into the equation. We also set n = 0.4 because this value 
seems to describe correctly the shear-thinning behaviour 
of our emulsion (see Figs. |l^ and |2(](b)). With only two 
free parameters <j and A, we were then able to fit the 
velocity profiles very closely so that the fits were undis- 
tinguishable from those obtained previously with Eq. ([is]) . 
However, once again, we could not find a single set of pa- 
rameters to fit all the velocity profiles. A and a had to 
be allowed to vary significantly with a so that the whole 
range of investigated shear rates could be accounted for. 

Due to the restricted amount of data and to the ex- 
perimental uncertainty, it was not reasonable to allow 
one more free parameter in the fits, namely to allow n 
to vary. In that case, almost any value of n between 0.1 
and 0.7 could be obtained by tuning A and a . Indeed, 
in the frawework of Eqs. ( |T§| ) and ([l9]), the curvature of 
v(x) is not only controlled by n but is also very sensitive 
to cr . However, if we restrict the analysis to the range 
0.2 < 7 < 5 s _1 where D(<r 1 ) appears as very linear on 
Fig. |l](b), we showed that n = 0.4 and <r = 90 ± 10 Pa 
were reasonable parameters for this range of applied shear 
rates. Note that this range of shear rates corresponds to 
the domain where the local flow curve of Fig. |l9| seems the 
most continuous. 



5.3.3 Existence of a global flow curve 



5.3.2 Influence of a yield stress 

In the case of the concentrated emulsion, it is very tempt- 
ing to interpret the particular value <j a of the stress where 
D(a ) — as a "yield stress". However, if we assume the 
following local behaviour based on the Hcrschcl-Bulklcy 
equation: 



cr(x) = a Q + Aj(x) n for a > <j a , 



(18) 



the above analysis of the velocity profiles is no longer valid. 
Indeed, if Eq. ( |l8| ) holds and if n ^ 1, the velocity v(x) is 
no longer easily integrable from Eqs. (||) and (||). Instead, 
for n^l, one gets: 



v(x) v 1 
R^+x~ R x 



Rx+x 



It 1 



dr 
r 



A 



l/n 



(19) 



Note that we can rule out the case n = 1 for two 
reasons: (i) although they do not allow for a reliable mea- 
surement of cr ; rhcological data indicate a power-law be- 
haviour with n w 0.4 at large 7 (see Fig. |l^), and (ii) n = 1 
leads to velocity profiles whose curvatures varie strongly 
from a logarithmic shape close to a to an almost linear 



The above problems encountered during the analysis of 
the individual velocity profiles raise the question of the 
existence of a global flow curve for the concentrated sys- 
tem. Indeed, it seems that simple approaches based on 
Eqs. (0) or (|l|) fail to describe our experimental data 
for the 75 % emulsion with a single set of parameters. 
The fact that the curvature of the individual profiles may 
be nicely fitted in terms of one of these two equations 
means that those flow behaviours may be representative 
of the fluid at a local level for a given shear rate. How- 
ever, except on a small range of 7, this kind of simple flow 
behaviour does not hold globally: the parameters have to 
be adjusted from one shear rate to another. This result is 
in agreement with measurements by Mason et al. (|] who 
report that for emulsions of droplet size 0.25 (im and con- 
centrated above 65 %, (7(7) deviates from a power law at 
high 7. 

Thus, our results would lead to a general picture of a 
sheared concentrated emulsion as a material whose flow 
behaviour changes with the imposed stress due to subtle 
changes in its micro-structure: 



a{x) =<7 a (a)+A{a)^t(x) n ^ 



(20) 



Indeed, even though we checked that the mean radius of 
the oil droplets did not change during our experiments, 
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the parameters controlling the flow behaviour such as the 
shear-thinning exponent n may depend on the compres- 
sion state of the droplets or on their precise shape or even 
on the sample history. 



6 Discussion and conclusions 

In this paper, we presented local velocity measurements in 
a dilute and in a concentrated emulsion using heterodyne 
DLS. In both systems, we found significant slippage. We 
have shown that the corrections for wall slip usually intro- 
duced in rheology are valid for the concentrated emulsion. 
However, in the dilute emulsion, incrtial effects may come 
into play and complicate the analysis so that direct esti- 
mate of slip velocities are required. 

Once wall slip is measured, we can access interesting 
local rheological information. We found that the dilute 
emulsion has a Newtonian behaviour and that the con- 
centrated system is shear-thinning with an exponent close 
to 0.4 that accounts well for the curvature of the velocity 
profiles. It also seems that the behaviour of the concen- 
trated emulsion cannot be described by a single flow curve 
over the whole range of shear rates that we investigated. 
This raises the question of the definition of a global flow 
curve. 

We have also shown that for the smallest shear rates 
at which the concentrated system flows (0.4 < 7 < 1 s _1 ), 
the flow behaviour is compatible with a yield-stress fluid 
with <t w 90 Pa and n ps 0.4. In any case, our results 
clearly rule out a simple power-law behaviour with a single 
prcfactor A (see Fig. 2l](b)). For 7 < 0.2 s" 1 , the emulsion 
does not flow and behaves like a solid. This shows the 
existence of two different regimes: one where the system 
is solid-like and one where the system flows like a yield- 
stress fluid. Up to the precision of the present data, the 
two shear rates at which the system is solid-like seem to 
correspond to shear stresses slightly above a Q (a = 100 
and 130 Pa, see Fig. |l^). This raises important questions 
about the behaviour of a concentrated emulsion near the 
transition between the two regimes. 

Indeed, recent experiments [| 18 and numerical simu- 
lations Q on soft disordered materials have focussed on 



the local flow properties of yield-stress systems like clay 
suspensions, industrial gels, emulsions and foams. Such 
systems were shown to present inhomogeneous flow pro- 
files when sheared close to their yield point: the sample 
separates between a fluid-like phase near the rotor and a 
solid-like phase near the fixed wall of a Couette cell. In 
some cases, the velocity profiles in the fluid phase show a 
strong curvature very similar to that observed on our mea- 
surements in the concentrated emulsion Q . These features 
seem quite general since they have also been observed in 
colloidal suspensions and granular pastes @,|l) . 

Since the first models of soft glassy materials under 
shear |2^,^3| and triggered by recent experimental devel- 
opments on aging and rejuvenation p4| . p5 26 , intense the- 
oretical effort is currently being spent on the understand- 
ing of the yield stress phenomenon [p7| , p8l|29| and on the 
possible existence of two separate branches in the flow 



behaviour |^,|3l|]. In particular, hysteretic effects are ex- 
pected when the shear stress is controlled and the presence 
of shear bands is predicted when the shear rate is imposed 
& 

So far, we have not been able to show that inhomo- 
geneous flows and/or hysteretic effects take place in the 
75 % emulsion when sheared in our small gap geometry 
(e = 3 mm for i?j = 22 mm). Direct visualization of the 
emulsion surface as in Refs. J2|,[l0| is made impossible by 
the use of a thermostated cover that prevents evaporation 
of our samples. Moreover, let us recall that our velocity 
profiles are obtained in typically 30 minutes. The flow pro- 
files of Fig. |^ showing solid-body rotation may thus miss 
transient inhomogeneous behaviours that could occur on 
shorter time scales such as fracture or stick-slip phenom- 
ena. 

However, none of the velocity profiles presented here 
seems consistent with a stationary banded flow as de- 
scribed in Refs. j|,[l8). More experiments around the yield 
point will be needed to conclude on this issue. In particu- 
lar, in the framework of the recent theoretical approaches 
of inhomogeneous flows, the existence of two points above 
the yield stress where the system behaves like a solid body 
should lead to inhomogeneous and possibly non-stationary 
flows in this region. Future experiments will focus on this 
specific point. We also intend to use high-frequency ultra- 
sound to measure velocity profiles with an increased tem- 
poral resolution (0.01-0.1 s) in order to capture transient 
and possibly non-stationary flows. 

The authors wish to thank the "Cellule Instrumentation" at 
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ette cells used in this study. We are also very indebted to 
B. Pouligny for his decisive participation to the optical setup 
and for helpful physical discussions. We thank A. Ajdari, C. Bar- 
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